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Abstract 



>• » 

q ' Approximate solutions representing the gravitational-electrostatic balance 

^ | ■ of two arbitrary point sources in general relativity have led to contradictory 

. arguments in the literature with respect to the condition of balance. Up to the 

present time, the only known exact solutions which can be interpreted as the 
\ non-linear superposition of two spherically symmetric (Reissner-Nordstrom) 

bodies without an intervening strut has been for critically charged masses, 
Mf = Qf . In the present paper, an exact electrostatic solution of the Einstein- 
Maxwell equations representing the exterior field of two arbitrary charged 
Reissner-Nordstrom bodies in equilibrium is studied. The invariant physical 
\Q \ charge for each source is found by direct integration of Maxwell's equations. 

The physical mass for each source is invariantly defined in a manner similar 
to which the charge was found. It is shown through numerical methods that 

i ■ balance without tension or strut can occur for non-critically charged bodies. 

?— i i 

(3JT), It is demonstrated that other authors have not identified the correct physical 

parameters for the mass and charge of the sources. Further properties of 
the solution, including the multipole structure and comparison with other 
jJJ ■ parameterizations, are examined. 

04.20.Jb, 04.40.Nr 
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I. INTRODUCTION 



In a recent paper by Bonnor M, the equilibrium conditions for a charged test particle 
in the field of a spherically symmetric charged mass (Reissner-Nordstrom solution) were 
investigated. He found that the classical condition for equilibrium 

M 1 M 2 = Q 1 Q 2 (1.1) 

for which the separation between the particles is arbitrary, was neither necessary nor suf- 
ficient for electrostatic balance of two spherical masses. This is in conflict with the earlier 
results of Barker and O'Connell and Ohta and Kimura [[| who used different approxima- 
tion methods. Barker and O'Connell claimed that in the post-Newtonian approximation, 
the equation 

(M 1 Q 2 -M 2 Q 1 )(Q 1 -Q 2 ) = (1.2) 



had to be satisfied in addition to (|1.1|) . Ohta and Kimura claimed that in the post-post- 
Newtonian approximation, the necessary and sufficient condition for balance is that each 
mass be "critically" charged, 

Mi=\Qi\ i = l,2 (1.3) 

and balance can occur for arbitrary separation of the sources. Up to the present time, the 
problem of gravitational-electrostatic balance of two spherical bodies in general relativity 
without an intervening Weyl line singularity (strut or tension) has been solved exactly only 
for critically charged masses A balance solution was originally thought to have been 

found |7| within the Herlt class for both sources having Mi > \Qi\, but it was subsequently 
shown that the intervening line singularity could not be removed ||. Kramer || presented 
an exact solution for the electrostatic counterpart of the double Kerr-NUT solution with 
zero spin parameter. He found that condition ( |1 . 1| ) holds for electrostatic balance. However 
he stated that his solution cannot be interpreted as the non-linear superposition of two 
Reissner-Nordstrom solutions and thus the masses are not spherically symmetric. 

In the present paper, an exact electrostatic solution of the Einstein-Maxwell equations 
representing the exterior field of two arbitrary charged non-linearly superposed Reissner- 
Nordstrom sources in equilibrium is given. It is obtained with the aid of Sibgatullin's flO| 



method for constructing the complex Ernst potentials [TT|. It is mathematically equivalent 
to the solutions of Manko et al JT^] and Chamorro et al |13|], henceforth referred to as 
papers I and papers II respectively, (with their spin parameters set to zero) and they are 



all special cases of the general mathematical solution given by Ernst It is of primary 
importance that the parameters in the solution be related to a physical set of parameters in 
order for any subsequent analysis of the solution to have any significant physical meaning. 
For a physical set of parameters, one would prefer to use the individual masses and charges 
of each source and the distance between the sources. The invariant charge enclosed by a 
spacelike hypersurface can be found by the direct integration of Maxwell's equations. For 
space-times with a timelike Killing vector, a conserved quantity which can be interpreted 
as the contribution to the total mass from each body can be invariantly defined in analogy 
with the charge (See, for example, Refs. |9|,|l5| and ||16|| ). This paper follows Kramer || for 
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the definition of the individual mass of each body. In Section ||, the integrals of charge and 
mass are given and they are applied to the Weyl-class solution for two Reissner-Nordstrom 
bodies in Section PTT[ Section [IV] presents the solution for a parameterization of the non- 
Weyl-class double Reissner-Nordstrom solution based on the Weyl-class parameterization. 
These are then compared to the parameterizations proposed in papers I and II. It is shown 
that the parameterizations employed in papers I and II do not represent the physical masses 
or charges of the individual sources even in the Weyl-class limit (except for the special case 
of identical bodies in paper I). Due to the complexities of the parameterization, a rendering 
of the solution in terms of the individual masses and charges as given in Section [IJ has not 
yet been accomplished. However, numerical analysis of the physical masses and charges is 
possible for a given set of parameters. In Section [V], balance without a strut or tension for 
numerical values of the physical mass and charge is examined. It is found that there are 
balance conditions for which neither body is critically charged and the Newtonian balance 
condition does not hold. This is in accordance with Bonnor's [|I| test particle analysis. The 
dependence of the balance condition on the separation of the bodies is not yet known. A 
discussion of the results and conclusions are given in Sections [VT| and |V11| . 



II. MASS AND CHARGE 

For a static axially symmetric space-time, the mass Mj and charge Qi of a source inside 
a closed 2-surface Oi are given by the integrals |T7|j9|| 

Mt = ~lK*J=gdr* (2-1) 

Qi = ~fF*V=9<fo (2-2) 

where 

R ab = ^a;b + ^pab _ ^.3) 

The timelike Killing vector is £ a , F ab is the electromagnetic field tensor, $ is the electrostatic 
potential, g is the determinant of the metric and df* b is the dual to the surface element 2-form 

df a \ 

df* ab = \e abcd df cd (2.4) 

(here e abc d is the flat space Levi-Civita permutation symbol). The above integral conserva- 
tion laws follow from the local conservation laws 

the first, following from the conservation of charge and the second from the existence of the 
timelike Killing vector £ a and the restriction to a static axially symmetric space-time metric. 
Since the Einstein-Maxwell equations also imply 

F ab . b = K ab b = 0, (2.6) 

in a source free region, any deformation of the surface <7j in the electrovacuum region outside 
the sources does not change the values of the integrals Mj and Qi . 
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III. THE WEYL-CLASS TWO BODY SOLUTION 



To investigate the structure of space-times with two sources, the Weyl-class double 
Reissner-Nordstrom solution provides a suitable yet mathematically uncumbersome frame- 
work from which to proceed. The solution is easily found through the method presented in 
Ref. The metric for a static axially symmetric space-time can be written in the canonical 
form 

ds 2 = e w dt 2 - e v ~ w (dp 2 + dz 2 ) - p 2 e~ w d<p 2 , (3.1) 

where w and v are functions of the cylindrical coordinates p and z. The Weyl-class solutions 
are characterized by the metric function w being a function of the electrostatic potential, 
i.e. w = u>($) so that the gravitational and electrostatic equipotential surfaces overlap. For 
asymptotically flat boundary conditions, the unique functional relationship between e w and 
$ is 



e w = l-2— $ + $ 2 , (3.2) 
q T 

where $ is the electrostatic potential and m T and q T are the total mass and charge re- 
spectively. The solution representing two "undercharged" (Mj > \Qi\) Reissner-Nordstrom 
bodies (or "black holes") is given by 

* = 4^1' (3 - 3) 



where 



= ( R 1 + R 2 -2l 1 \ ( R 3 + R 4 -2l 2 \ 
J \R 1 + R 2 +2l 1 [ R 3 + R 4 + 2l 2 ' l ' j 



R\ = (z - d - 2k) 2 + p 2 , 
R\ = {z- d) 2 + p 2 , 

(z + d + 2l 2 ) 2 + p 2 . 



R 
R 2 



(3.5) 
(3.6) 
(3.7) 
(3.8) 



The constant parameters 2d and 2li,2l 2 are the coordinate distance between the horizons 
and the "lengths" of the horizons (Weyl "rods") respectively (see Fig. [1]). The parameter a 
is defined through the equation 



2m-, 



(3.9) 



The metric function e w is found through Eq. 



(Ri + R 2 f - Al 2 (R3 + R4) 



til 



4:R\R 2 



The metric function e v is 



((h + l 2 + d)R 1 + (l 2 + d)R 2 - hRj) d 
((Zi + d) Rt + R 2 d - hR 3 ) (Z 2 + d) 



Choosing the surface o\ to encompass body 1 and the surface a 2 to encompass body 2 of 
Fig. |l], the mass and charge integrals of Eqs. (|2.1| ) and (|2~2l) yield 



1 + a 2 , ^ 2a , 

M 1 = - -h Q 1 = -h 

1 — a 1 1 — a 1 

(3.10) 

1 + a 2 , ^ 2a , 

M 2 = -l 2 Q 2 = -l 2 

1 — a 1 1 — a 1 

The above form of the individual mass and charge for each Reissner-Nordstrom body is 
similar to the form proposed in Ref. 0] for the mass and charge decomposition of two charged 
Curzon particles. It was stated in Ref. |7j that the conjectured charge decomposition for 
both the double Reissner-Nordstrom and Curzon cases were verified by direct calculation 
through Eq. ( |2.2|) . It is straightforward to verify that Eq. ( |2.1| ) yields the conjectured 
mass decomposition for the double Curzon solution. Because of the functional relationship 
between the gravitational potential and the electrostatic potential, not all of the parameters 
Mi, M 2 , Qi, Q 2 are independent. Thus the Weyl-class is also characterized by the constraint 

M 1 Q 2 = M 2 Q 1 . (3.11) 

Removal of the line singularity between the bodies yields Eq. as an additional condition 
on the parameters. As a result, the parameters also satisfy Eq. ( |1.3|) . Thus equilibrium 
without a strut or tension occurs for "critically" charged sources and this balance is found 
to be independent of the separation distance M. 



IV. NON-WEYL PARAMETERIZATIONS 

Generalizing the Weyl-class double Reissner-Nordstrom solution to the case in which 
the gravitational and electrostatic equipotential surfaces no longer overlap has usually been 
attempted through the means of generating techniques (see, for example, Refs. ||[7| and ||). 
In these techniques, new solutions are generated from old ones rather than by solving the 
equations directly. Recently, considerable interest has focused upon a method []T0|| which 



constructs the Ernst potentials |IT] from initial data on the symmetry axis. The complex 
Ernst potentials £(p, z) and ^(p, z) of all stationary axisymmetric electrovacuum space-times 
with axis data of the form 

where 

U = z 2 + U lZ + U 2 (4.2) 

V — V\z + V 2 (4.3) 

W = Wiz + W 2 (4.4) 
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and Ui,U2,Vi,V2, Wi, W 2 are complex constants, have been found []I4| . However, a mathe- 
matical solution to the Einstein-Maxwell field equations does not imply a well understood 
physical interpretation of the solution. Sibgatullin's method of constructing the Ernst po- 
tentials aids in obtaining the physically meaningful parameterization which is sought for the 
two-body case in question. 

In Sibgatullin's method, it is required that the Ernst potentials along the z-axis be 



specified. Our choice was 



S(p = 0,z) = e(z) 



1 - 



2(mi(z + z 2 ) + m 2 {z + Zi)) 
[z + Z\ + mi) (z + z 2 + m 2 ) - q\q 2 ' 



,w n \ 771/ \ qi(z + z 2 ) + q 2 (z + Zl ) 

W(p = 0, z) = F[z) = — r . (4.5) 

[z + Zi + mi)(z + z 2 + m 2 ) - q x q 2 

It has the form of the Weyl-class double Reissner-Nordstrom axis data. If the additional 
Weyl-class constraint 

niiq 2 - m 2 qi = (4.6) 

is placed on the functions e(z) and F(z), then Sibgatullin's method yields the Weyl-class 
double Reissner-Nordstrom solution (in an alternate form to Ref. Q) and the parameters 
mi,m 2 ,qi,q 2 are the physical masses and charges as defined by Eqs. (pTTH2.2|) (i.e., Mi = 
rniiQi = qi,M 2 = m 2 ,Q 2 = q 2 ). For the solution of two Weyl-class Reissner-Nordstrom 
black holes (given in Section p|) , Fig. [3] shows the coordinate positions of the centers of 
the "rods" as d + l\ for body 1 and — d — l 2 for body 2. The parameters z%, z 2 identify the 
negative of the coordinate positions of the centers of the "rods", i.e., 

Z\ = — d — l\ , 
z 2 = d + l 2 . 

If condition ( |4.6| ) is not imposed, w 7^ w(Q), i.e., the gravitational and electrostatic equipo- 
tential surfaces no longer overlap. In Section it will be shown that the parameters 
mi,m 2 ,qi, q 2 then no longer carry the suggested physical meaning and the parameters z%, z 2 
no longer coincide with the centers of the "rods" when the Weyl-class constraint (|4.6|) is not 
imposed. 

The full Ernst potentials £(p, z) and \l/(p, z) for the axis data of Eq. (f!5|), expressed in 
terms of the cylindrical coordinates (p, z), are found to be (the details of the method can 
be found in Refs. f^OlOBll and in the review article IE21) 



£ = ^ ¥ = (47) 



where 

4 44 



A^^aij^rj, B = J2 b i r i^ C = J2°i r ^ r n = ^ p 2 + (z - a n ) 2 , (n = 1 -> 4). 

i<j i=l i=l 

The constants a n in Eq. (|4.?] ) are the roots of the equation 
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e(z) + [F{z)f = 



(4.8) 



and can only be real or complex conjugate pairs. The remaining constants a^, hi and q are 
defined as follows: 



SkV k SiVi 
t k U k tiUi 



(i <j; k <l- k,ly£ i, k = 1 -> 3; j,l = 2 -> 4); 



"1) Si^iySi 



„2 f 2 



2*2 



s 2 t 2 



(k < I <m; k,l,m ^ i; i = 1 — > 4; fc = f , 2; / = 2, 3; m = 3,4); 

Cj = i-iy+hMsi - u)(K 3 Gi + k 4 ^), 



Sktl 



sit\ s m t m 



s m v 1 

t k U k tiUi t m u m 



Hi 



s\H k sfti S^tm 
S k V k SiVi s m v m 
t k U k tiUi t m u m 



(4.9) 



(k <l < m; k,l,m^i; i — 1 — > 4; fc = 1,2; / = 2,3; m = 3,4); 



Si = Pi - Q!j, U = f3 2 - Cti 



_ m 1 z 2 + m 2 z 1 + {mi + m 2 )(3i T , _ m x z 2 + m 2 z l + (m 1 + m 2 )(3 2 

K l = n n > K 2 = 



Pi -Pi 



- <li z 2 + 92*1 + (qi + q 2 )Pi _ q\z 2 + q 2 z x + {q x + q 2 )p 2 

K 3 = n n > K 4 = 



Pl-P2 



P2-Pl 



Pi = (^1 + m 1 + z 2 + m 2 - \J (z! - z 2 + m 1 - m 2 f + Aq^^j , 



P2 = (zi + m 1 + z 2 + m 2 + \J (z x - z 2 + m 1 - m 2 f + Aq x q 2 ^j , 

where all of the subsequent quantities introduced are constants ultimately defined in terms of 
rrii, q iy z iy % — 1,2, which specify the character and locations of the sources in the Weyl-class 
limit only. 
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The expressions for 8 and \I/ are in Kinnersley's [23 form and this permits one to write 
the corresponding metric functions as 



AA-BB + CC AA-BB + CC 

" (A + B)(A + BY K ori r 2 r 3n ' 1 ' 1 

where 

\Kj J \i<j J 

and a bar denotes complex conjugation. For a static metric, the electrostatic potential $ is 
equal to the Ernst potential ^ and this completes the solution. 

With the knowledge of the full Ernst potentials and the metric functions, the next 
step would be to evaluate the true mass and charge integrals in terms of the parameters 
mi, m,2, qi, qi, %i, It is to be stressed that outside of the Weyl-class, these parameters no 
longer carry the suggested physical meaning. For the metric (|3.1|) , the integrals (|2.1| ) and 
can be written as relations in flat 3-space (i = 1,2): 

Mi = — £ e~ w £ a n a dA (4.12) 
8n Jui 



Qi = --?- <fe~ w $, a n a dA, (4.13) 

where n a (a runs from 1 to 3) is the unit vector orthogonal to the surface and dA denotes 
the invariant (flat) surface element (see also Ref. || and references therein). 

We can extend the Weyl-class definitions of the coordinate positions of the bodies to 
the non- Weyl-class solution. There are three distinct types of sources of interest. They 
are characterized by the transition between a source with an event horizon to one without 
an event horizon. As mentioned previously, the constants a n , n = 1 — > 4 in Eq. ( |4.7| ) 
are either real or complex conjugate pairs. By definition, we choose a\ > a 2 > Q3 > 
«4. A Reissner-Nordstrom "black hole" is characterized by real pairs of a n . Fig. [I] shows 
that, in the Weyl canonical coordinate system, a n indicates the end points of a "Weyl 



rod", which itself is the event horizon surface. A "superextreme" object |12j or "naked 
singularity" is characterized by a complex conjugate pair of a n . Body 2 of Fig. || illustrates 
the manifestation of a "superextreme" body in the space-time. An "extreme" object, for 
example, would be characterized by real a n for which a\ = a 2 - Therefore we have the 
following definitions for the coordinate positions of the sources: 

i) For a Reissner-Nordstrom "black hole", we define — Z\ to be the coordinate position of 
the center of the "Weyl rod". For example, the coordinate position of body 1 of 

Fig. ® is 



s 



ii) For a "superextreme" object, we define — Zi to be the coordinate position of the real 
part of a n . For example, body 2 of Fig. (H) is a "superextreme" object. Therefore its 
coordinate position is 

— Z 2 = Re(a 3 ) = Re(a4). 

(One could consider the imaginary part of a n as the end points of a "complex Weyl 
rod" with the coordinate position of this "complex rod" being defined as its intersection 
with the real axis (z-axis).) 

Hi) For an "extreme object", we define — Zi to be the coordinate position of the point 
locating the zero "length" Weyl "rod". For example, if body 1 was an "extreme" 
object, then a.\ = a 2 and —Z\ = c^. 

We also define 

Re(a 2 ) > Re(a 3 ) (4.14) 

as the condition for having two separated bodies irrespective of the type of object. 

With the above integrals and the coordinate positions as defined above evaluated in terms 
of mi, m.2, qi, q 2 , z ii z 2, it would then be possible, in principle, to invert these equations 
and hence write the solution fl4.7[ - |4-9|) in terms of the true physical parameters Mi, Qi 



and the coordinate positions Zi, i — 1,2. Ideally, the coordinate positions of the sources 
should be replaced with the proper separation of the sources. The complexity of the above 
Ernst potentials makes the analytic evaluation of the integrals ( 4.12 - pTT5 ) and the proper 



separation difficult. As a consequence this goal has not yet been achieved. However, it is 
possible to numerically integrate Eqs. (PTT|-|2.2D for a given set {mi,m 2 , q%, q 2 , Z\, £ 2 }. This 



will prove to be useful in studying balance conditions without a strut in Section [V[ 

Although the numerical evaluation of the physical mass and charge can be achieved from 
the parameterizations of paper I or paper II, it was hoped that the parameterization proposed 
in this paper, based on the Weyl-class solution, would facilitate the analytic evaluation of 
the integrals. It is not difficult to show that the parameterizations in papers I or II do 
not correctly identify the individual masses and charges of each source. We stated earlier 
that our parameterization {mi, 777,2, qi, q%, Z\, z 2 } only represents the physical masses and 
charges and coordinate positions of each source when the Weyl-class condition (Eq. ( |3 . 2|) 
or flop) is imposed (i.e., {M x = m 1 , M 2 = m 2 ,Q\ = qi,Q 2 = <?2,-Zi = z 1 ,Z 2 = z 2 }) ■ 
We can best demonstrate the problems with the parameterizations of papers I and II by 
comparing the representation of a properly parameterized Weyl-class solution with each of 
the other parameterizations. Let the set {mi, m 2 , qi, q 2 , zi, z 2 } represent the physical Weyl- 
class parameters under the condition mig 2 = m 2 q\. Then the relationships between the 
three parameterizations is found by solving the set of equations (setting the spin parameters 
found in papers I and II to zero) 
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Weyl-class Paper I Paper II 

mi + m,2 = rhi + "^2 = Th\ + m 2 

Q1 + Q2 = <?i + <?2 = qi + Q2 

zi + z 2 = Z1 + Z2 = Zi + z 2 (4-15) 

rriiZ2 + ^2^1 = fhiz 2 + m^i = rhiZ 2 + m 2 Zi + 1ih\m2 

q x Z 2 + q 2 Z\ = q\z 2 + q 2 Zi = q\Z 2 + q 2 Zi + qim 2 + q 2 rhi 

Z\Z2 + m 1 m 2 - qiq 2 = Ziz 2 + mira 2 = z x z 2 - mim 2 . 

The tilded and careted parameters are the parameterizations of papers I and II respec- 
tively. Table [I] summarizes the results of solving the system ( |4.15| ) given the values shown in 
column 1. The solution represents two Weyl-class Reissner-Nordstrom "critically charged" 
bodies without an intervening line singularity. It is clear that none of the parameter values 
in the latter two columns match the physical Weyl-class values. In fact one has to assign 
negative values to m,2, q_2 in order to obtain a positive physical mass and charge for source 2. 
Thus, apart from one special case, neither the paper I nor the paper II parameterizations can 
be interpreted in as the invariant physical parameters. The only exception is for identical 
bodies (with or without a line singularity) in the parameterization of paper I. In this very 
special case of the Weyl-class, the parameters rhi = m 2 , qi — q\ are the physical masses 
and charges. However, z\ and z 2 do not identify the coordinate positions of the bodies as 
defined earlier. The paper II parameterization is not physical even for identical bodies. 

It is the demand for the inclusion of the Weyl-class solution in Ref. j|] which led to 
our form of e(z) and F(z). It should be emphasized that our parameterization contains 
as a special case, the simplest clearly individually spherical two-body balance solution of 



two critically charged bodies. This can be best illustrated by examining the Simon p4| . [25 
relativistic multipole moments of each parameterization. The first five Simon relativistic 
mass and charge multipole moments for our parameterization are 

Mo = mi + m 2 , 

Mi = - m 1 z 1 - m 2 z 2 , 

M 2 = mizf + m 2 z\ - (m 1 m 2 - qiq 2 ) {mi + m 2 ) , 

M3 = - m 1 zf - m 2 zl + (7711TO2 - gi<?2) (2mizi + 2m 2 z 2 + zim 2 + z 2 m 1 ) , 
Ma = mizf + m 2 z\ - (mim 2 - qiq 2 ) (mi + m 2 ) (qiq 2 - m 1 m 2 ) 

+ 2 ( mi zf + m 2 z\) + (mi + m 2 ) {z 1 + z 2 f ( 4 - 16 ) 

+ i (mi + m 2 ) ((gi + q 2 f - (mi + m 2 ) 2 ) 

~ 2T0 ( Zl ~ 22 ) 16 ( Zl ~ z ^> ( mi + m2 ) ( mi< ? 2 _ m 2gi) 2 

+ z\ (30mi (mim 2 + m\ - q 2 ) - 3gi (3m 2 gi + 7q 2 mi)) 
- z 2 (30m 2 (mim 2 +m\ — q\) - 3q 2 (3miq 2 + 7qim 2 )) 

and 
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(4.17) 



Qo = qi + q2, 
Qi = - q x zi - q 2 z 2 , 

0,2 = qizf + q 2 z\ - (m 1 m 2 - qiq 2 ) (qi + ^2) , 

Q3 = - qizf - q2z\ + (m 1 m 2 - qiq 2 ) (2q 1 z 1 + 2q 2 z 2 + z x q 2 + z 2 qi) 
Qa = q\z\ + q 2 z\ - (m 1 m 2 - qiq 2 ) (q x + q 2 ) (qiq 2 - m x m 2 ) 
+ 2 {q x z\ + q 2 zl) + (g x + q 2 ) \z x + z 2 f 

+ \ (qi + <?2) ((gi + q2) 2 - (mi + m 2 ) 2 ) 

~ 2l0 ( 2l ~ Z ^> [ 16 ( Zl ~~ ^2) + 92) ( mi ^ 2 ~~ m 2<?l) 2 

- (30g 2 (?i?2 - ^1^2 + ?i) - 3mi (13mi^ 2 - 3m 2 <2i)) 
+ z 2 (30gi (gi^ 2 - »im2 + qf) - 3m 2 (13m 2 gi - 3mig 2 )) 

respectively. In Newtonian physics, a system of two monopoles at positions zi, z 2 has 
multipole moments 

M n = m lZ ^ + m 2 z^, Qn = qiz? + q 2 z2. (4.18) 

It is interesting to observe that this is also the relativistic multipole structure for two Weyl- 
class critically charged bodies, at least up to Ma, Qa- There is an inherent asphericity 
imposed upon each, since the two bodies are interacting in a line. For non-linearly interacting 
sources in a line, one would not expect to realize perfect sphericity of the individual sources. 
(It is yet to be explained why the sphericity is maintained in the Weyl-class, at least up to 
M4, Qa-) Once the solution is written analytically in terms of the physically meaningful 
constants Mi,Qi and the coordinate positions Z i} i — 1,2, one will be able to examine the 
general multipole structure of non-linearly interacting spherical bodies. 

For comparison, the first four Simon relativistic mass and charge multipole moments for 
the parameterization of paper I (with their spin parameters Oj = 0, i — 1, 2) are 

Mo — rhi + m 2 , 

Mi = - rhxlx - rh 2 z 2 , (419) 
■■ m\z\ + m 2 z\ - m x m 2 (mi + rh 2 ) , 



and 



M 2 

M3 = — rh\z\ — m 2 z\ + m\m 2 (2m 1 zi + 2m 2 z 2 + z\m 2 + z 2 rh\) 



Qo = qi + <?2, 

Qi = - q\z\ - q 2 z 2 , 

Q2 = q\z\ + q 2 z\ - m x m 2 (q x + q 2 ) , 

Q3 = ~ qizf - q 2 z\ + m x m 2 {2q x z x + 2q 2 z 2 + z x q 2 + z 2 q x ) 



(4.20) 



The first four Simon relativistic mass and charge multipole moments for the parameterization 
of paper II (with their spin parameters a, = 0, % = 1, 2) are 

Mo = m 1 + rh 2 , 

M\ = — rh\Z\ — rh 2 z 2 + 2m 1 m 2 , 

M 2 = rh\z\ + m 2 z\ + rhxm 2 {fh\ + m 2 — 2i\ — 2z 2 ) , (4-21) 
M3 = — ih\z\ — m 2 z\ 

+ m\m 2 {2rh 1 m 2 + 2z±z 2 + 2z\ + 2z\ — rhiz 2 — rh 2 zi — 2rh\Zi — 2m 2 z 2 ) 
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and 



Qo = gi + 92, 

Qi = - qyz\ - q 2 z 2 + rh x q 2 + m 2 q\, 

Q 2 = q\z\ + q 2 z\ + mim 2 {q\ + q 2 ) - {qifn 2 + fethi) (zi + z 2 ) , (4.22) 
Q 3 = - (ftif - g 2 5| - mim 2 (2^1^ + 2g 2 i 2 + ^ig 2 + £ 2 <7i) 
+ (qim 2 + q 2 rhi) (m 1 m 2 + z\Z 2 + z\ + i 2 ) • 

If the above parameterizations did represent the physical mass and charge, it is evident that 
the multipole structure would not be that of Newtonian spherical bodies even for critically 
charged bodies. As stated earlier, it should be noted that in the parameterization of paper I, 
it can be shown that only in the case of identical bodies, the parameters rh\ = rh 2 , qi = q 2 
are the physical mass and charge. However, in this case the multipoles still do not have the 
form of Eq. ( f4.18p since the parameters z\ and z 2 do not identify the positions of the bodies 
as defined earlier. A simple transformation would correct the multipoles in this case. 



V. THE EQUILIBRIUM CONDITION 

In order for the space-time to be regular on the z-axis between the sources (removal of 
the Weyl line singularity or imposition of the condition for elementary flatness ||26||), it is 
required that the metric function 

v ( Zi p = 0) = (5.1) 

between the sources. If the origin of the coordinate system is located between the sources 
(i.e., Re(ct 2 ) > 0, Re(«3) < 0), then application of Eq. (|5.1|), after some simplification, yields 
the balance equation 

r , a 12 (013 + a u ) + a 12 (a 13 + a u ) n oX 
K = = 0. (5.2) 

Three cases were examined: i) Two Reissner-Nordstrom black holes, ii) Two Reissner- 
Nordstrom superextreme bodies and Hi) One black hole and one superextreme body. 

The procedure for testing for equilibrium without an intervening strut or tension will be 
as follows: 

1. Assign numerical values to five of the six parameters from the unphysical set 
{m 1 ,m 2l q 1 ,q 2 ,z 1 ,z 2 }. 



2. Solve Eq. ( |5.2|) for the unknown variable. 



3. If a real root of Eq. ( |5.2j ) exists, then evaluate Eqs. ( |2.1|) and ( |2.2| ) to determine the 
physical mass and charge parameters. 

The results for each of the three cases are as follows: 
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A. Two Reissner-Nordstrom Black Holes 



Numerous sets of the parameters {mi, m 2 , q%, q 2 , Zi, z 2 }, such that the constants a n , n = 
1 — >• 4 are real, were investigated. No roots were found of Eq. (|5.2|) . For example, choosing 
mi = 9.0, qi = 3.0, Z\ = —15.0, m 2 = 8.0, q 2 = 2.0, no balance for < z 2 < 10 10 was 
found. These findings are consistent with other results PfljlB]] that two Reissner-Nordstrom 
black holes cannot be found in equilibrium without an intervening strut or tension. 



B. Two Reissner-Nordstrom Superextreme Bodies 

Numerous sets of the parameters {mi, m 2 , qi, q 2 , z±, z 2 }, such that the constants a n , n = 
1 — ► 4 are complex conjugate pairs, were investigated. No roots were found of Eq. (|5.2|). For 
example, in choosing m x = 3.0, q x = 9.0, z\ = —15.0, m 2 = 2.0, q 2 = 8.0, no balance for 
< z 2 < 10 10 was found. These findings suggest that two Reissner-Nordstrom superextreme 
bodies cannot be found in equilibrium without a strut or tension. 



C. One Black Hole and One Superextreme Body 



The following three different cases were found for which Eq. ( |5.2| ) has a real root. Each 
case has the configuration illustrated in Fig. [2|. 

Case A) For mi = 6.0, q\ = 2.0, z\ = —5.0, m 2 = —0.7, q 2 = 4.0, balance at 
approximately |^7j z 2 = 2.08 was found. The values of a n are a\ = 10.3, a 2 = 
1.74, a 3 = -3.11 + z4.30, a 4 = -3.11 - ?4.30. Using Eqs. (P^P^D , the physical 



masses and charges are M\ = 3.95, Q\ = —0.887, M 2 = 1.35, Q 2 = 6.89. 



Using the definitions of coordinate positions described in Section [TV], it was found 
that Z\ = —6.03 and Z 2 = 3.11. Thus balance has occurred for MiM 2 > QiQ 2 , 
QiQ 2 < at a coordinate separation of S = Z 2 — Z 1 = 9.13 . Note that the parameter 
m 2 is negative but both physical masses are positive. The parameterizations of papers I 
and II yield respectively 



Paper I 
mi = 4.96 
qi = 2.31 
m 2 = 0.34 
q 2 = 3.69 
z x = -6.60 
z 2 = 3.68 



Paper II 
rh\ = 4.36 
qi = -1.05 
m 2 = 0.94 
q 2 = 7.05 
z x = -6.00 
z 2 = 3.08 



which do not agree with the integrated values of Eqs. fl2.1|) and (|2.2|) . This demon- 
strates that in general none of the analytic parameterizations proposed, including our 
own, are suitable choices for the individual masses and charges of the sources. 

Case B) For mi = 9.0, qi = 3.0, z\ = —40.0, m 2 = 2.5, q 2 = 8.0, balance 
was found at approximately z 2 = 34.6. The values of a n are a% = 48.4, a 2 = 
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31.61, «3 = —34.62 + z7.65, «4 = —34.62 — z7.65. The physical masses and charges 
are Mi = 8.87, Qi = 2.00, M 2 = 2.63, Q 2 = 9.00. The coordinate positions are 
-Z x = 40.01, -Z 2 = -34.6 Thus balance has occurred for M 1 M 2 > QiQ 2l QiQ 2 > 
at a coordinate separation of S = 74.6 . 

Case C) For m x = 900.0, q l = 300.0, = -865.0, m 2 = 0.025, g 2 = 0.080, balance 
was found at approximately z 2 = 21.581. The values of a n are ct\ = 1713.5, a 2 = 
16.474, a 3 = -21.582 + i0.26226, a 4 = -21.582 - i0.26226. The physical masses 
and charges are M x = 899.71, Q x = 298.25, M 2 = 0.31897, Q 2 = 1.8254. The 
coordinate positions are —Z\ = 865.00, — Z 2 = —21.582 Thus balance has occurred 
for QiQ 2 > MiM 2 , QiQ 2 > at a coordinate separation of S = 886.58 . 

D. Comparison with Test Particle Analysis 

Bonnor's [|T| examination of a test particle in the field of a Reissner-Nordstrom source 
yielded a wide variety of balance conditions. The following cases for separation-independent 
equilibrium were examined (note: M, Q characterize the Reissner-Nordstrom space-time 
and m, q are the test body parameters): 

Case 1) For q = em, Q = t]M, e, i] = ±1, balance occurs if e = rj. 

Case 2) If m — \q\, \Q\ , or m ^ \q\, M = \Q\ , no equilibrium is possible. 

Case 3) If mM = qQ but m ^ \q\, then no equilibrium is possible. 

Since the exact solution under study contains the Weyl-class solution as a special case, we 
also find Bonnor's separation-independent equilibrium condition. Case 2 or 3 

cannot be tested readily by our numerical procedure. In order to do so, one would have 
to have the good fortune of correctly choosing the set {mi,m 2 ,qi,q 2 , z±, z 2 } such that the 
physical masses and charges satisfy the given conditions (i.e. Mi = \Q\\ etc.). Then, to test 
the dependence on separation, one would need to choose a new set of unphysical parameters 
such that the proper separation changes while the physical masses and charges remain the 
same. 

The following separation-dependent cases were also found in Ref. [l|]: 
Case 4) If \Q\ > M, mM = —qQ and m? ^ q 2 with qQ < 0, then an equilibrium exists at 

Q 2 

r = . 

2M 

Case 5) If \Q\ > M, \q\ < m, qQ < or 

Case 6) if \Q\ > M, \q\ < m, qQ > 0, qQ < mM then an equilibrium position exists at 

Q 2 (M [m 2 - q 2 ) + q^(m 2 - q 2 ) (Q 2 - M 2 )) 
r ~ m 2 M 2 - q 2 Q 2 
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Case 7) If \Q\ < M, \q\ > m, qQ > 0, qQ > mM then an equilibrium position exists at 

Q 2 (M (m 2 - q 2 ) - q^J{m 2 - q 2 ) (Q 2 - M 2 )) 
r = m 2 M 2 _ q 2Q2 • 

Thus we have found a direct correspondence between cases A-C of the exact solution and 
cases 5-7 of Bonnor's test particle analysis. The separation dependence of cases 4-7 cannot 
be studied in the exact solution using the present methods for the same reasons cases 2-3 
cannot be studied. Since the separation dependence cannot be tested using the present 
methods, there is little value in numerically calculating the proper separation of the sources 
in cases A-C. 

The physical parameters in case C could approximate a test body in a strong gravitational 
field. Using these values in case 7 and transforming from spherical coordinates to cylindrical 
coordinates for a single Reissner-Nordstrom body using the transformation (with 9 = 0) 

z = (r — M) cos 9, 

p = ^r 2 - 2Mr + Q 2 sin 6, (5.3) 

Bonnor's method yields a coordinate separation of <S = 1465.5. Since the separation of the 
bodies from these two methods are not consistent, it would appear that case C does not 
sufficiently approximate a test body. 



VI. DISCUSSION 

The essential departure in the present paper from previous work is the attempt to pa- 
rameterize the solution in terms of true physical constants of the space-time. For a static 
axially symmetric solution of the Einstein-Maxwell equations, the integrals of Eq. (|2.1| ) and 
( p.2| ) provide the invariant parameters required for meaningful analysis of the properties of 
the solution. 

There are three cases of the exact solution which have not been examined. They are 
an extreme body with respectively a Reissner-Nordstrom black hole, a superextreme body, 
and another extreme body for which the solution is not of the Weyl-class. Knowledge of the 
solution analytically in terms of the physical parameters is required to analyze these cases 
adequately. 

Ref. H defines the terms "undercharged", "overcharged" and "critically charged" as 
follows (i = 1,2): 

M 2 > Q 2 "undercharged" (6.1) 
M 2 < Q 2 "overcharged" (6.2) 
M 2 = Q 2 "critically charged" (6.3) 



For the Weyl-class, the "lengths" of the Weyl rods are § 2i< = 2JM 2 - Q 2 , i = 1, 2. If body 



1 is "critically charged" ||28|| , then a\ = a-i (= d) since 1% = (see Fig. |T]). This implies that 
the terminology "critically charged" body and "extreme" body may be used interchangeably 
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for Weyl-class solutions. If body 1 is "undercharged", ot\ (= d + 2li) and «2 (= d) are real 
quantities. Thus "undercharged body" and "black hole" are synonymous terms in the Weyl- 
class. Finally, if body 1 is "overcharged", «i (= d + li) and a 2 (= d + li) are complex 
conjugates. Thus the terms "overcharged" and "superextreme" are equivalent descriptions 
in the Weyl-class. Unlike the Weyl-class solutions where the "lengths" of the Weyl rods (real 
or complex) depend only upon the mass and charge of that source, it is strongly suggested 
from the analysis of Section [V] that for the general (non- Weyl-class) solution, the "lengths" of 
the rods also depend on the mass and charge of the other source and the distance separating 
the bodies as well. It would thus be possible to have a "critically" charged body (according to 
Eq. ( |6.3| )) for which the "rod" is either of non- zero "length" or "complex". This is important 
in terms of nomenclature for describing the physics of the space-time. Since the transition 
of a pair (eg. (ai,«2)) from real values to a complex conjugate pair in Sibgatullin's |1C 



method defines a differentiation of an object with a horizon to one without, it would seem 
that the appropriate description would be respectively, a black hole (horizon), "extreme" 
body (zero "length" Weyl rod) and "superextreme" body (no horizon or naked singularity) 
as described in paper I. The descriptions "under", "over" and "critically" charged body 
should be reserved for the relations Mf > Qf, Mf < Qf and Mf = Qf respectively between 
the individual masses and charges. This classification scheme would describe equilibrium 
conditions more precisely once all are identified. The appropriateness of such a scheme would 
become apparent when the analytic physical parameterization of the solution is known. 

Bonnor's |l| test particle analysis has been modified |30| in such a way that the equilib- 
rium conditions of a charged test particle in the field of a Kerr-Newman source can be stud- 
ied. The generalization of the mathematical solution to two spinning sources (Kerr-Newman 



sources) is already known |T^|T^|. One is able to invariantly define angular momentum for 
a stationary space-time in a manner similar to Eqs. ( |2.1| - |2.2| ) because of the presence of a 
spacelike Killing vector (rotational symmetry) (see Ref. |15| and references therein for def- 
initions of mass and angular momentum of stationary vacuum fields). It is unknown how 
the subsequent analysis of two identical spinning bodies in paper I based on the invariant 
definitions will affect their results, if at all. However, it is clear that the parameterization 
given is inadequate for the physical analysis of the general case (non-identical bodies). 



VII. CONCLUSIONS 

The solution derived in papers I, II and this paper is a generalization of the Weyl-class 
double Reissner-Nordstrom solution. However, the analytic parameterizations presented in 
papers I, II and this paper cannot in all cases be interpreted as the true physical constants of 
the spacetime. The invariant physical charge for each source is found by direct integration 
of Maxwell's equations. The physical mass is invariantly defined || in a manner similar to 
which the charge was found. Numerical methods were used to evaluate the invariant individ- 
ual masses and charges for the axially symmetric superposition of two Reissner-Nordstrom 
bodies. It was found that neither the Newtonian balance condition nor critically charged 
bodies are necessary for electrostatic equilibrium. The dependence of the balance condition 
on the separation of the bodies is not yet known due to the complexities involved in express- 
ing the solution analytically in terms of the true physical set of parameters. However, all the 
balance conditions found are consistent with Bonnor's test particle analysis. This suggests 
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that there exist equilibrium conditions which depend on the separation of the sources. The 
parameterization of this paper is manifestly physical in the Weyl-class limit. 

ACKNOWLEDGMENTS 

This research was supported, in part, by a grant from the Natural Sciences and Engineer- 
ing Research Council of Canada and a Natural Sciences and Engineering Research Council 
Postgraduate Scholarship (GPP). 



17 



REFERENCES 



W.B. Bonnor, Class. Quantum Grav. 10, 2077 (1993). 

B. M. Barker and R.F. O'Connell, Phys. Lett. A 61, 297 (1977). 
T. Kimura and T. Ohta, Phys. Lett. 63A, 193 (1977). 

F. I. Cooperstock and V. de la Cruz, Gen. Rel. Grav. 10, 681 (1979). 

A. Papapetrou, Proc. R. Irish Acad. 51, 191 (1947); S.M. Majumdar, Phys. Rev. 72, 
390 (1947). 

A. Tomimatsu, Prog. Theor. Phys. 71, 409 (1984). 

J. Carminati and F.I Cooperstock, Class. Quantum Grav. 8, L171 (1991); Gen. Rel. 
Grav. 24, 881 (1992). 

G. P. Perry, in Proceedings of the 5th Canadian Conference on General Relativity and 
Relativistic Astrophysics, R.B. Mann and R.G. McLenaghan eds., (World Scientific, 
Singapore, 1994) p. 449. 

D. Kramer, Class. Quantum Grav. 5, 1435 (1988). 

N.R. Sibgatullin, Oscillations and Waves in Strong Gravitational and Electromagnetic 
Fields (Nauka, Moskow, 1984), in Russian [English translation: Springer- Verlag, Berlin, 
1991]. 

F.J. Ernst, Phys. Rev. 168, 1415 (1968). 

V.S. Manko, J. Martin, and E. Ruiz, J. Math. Phys. 35, 6644 (1994). 

A. Chamorro, V.S. Manko and N.R. Sibgatullin, Lect. Notes Phys. 423, 119 (1993). 

F. J. Ernst, Phys. Rev. D 50, 4993; 6179 (1994). 

W. Dietz and C. Hoenselaers, Ann. Phys. (USA) 165, 319 (1985). 

H. Stephani, General Relativity, (Cambridge University Press, Cambridge, 1990), 2nd 
ed., p. 208 

Notations and conventions used are those of L. D. Landau and E. M. Lifshitz, The 
Classical Theory of Fields (Pergamon Press, Oxford, 1975), 4th revised English ed. 
H. Weyl, Ann. Phys. (Leipzig), 54, 117 (1917). 

G. P. Perry and F.I. Cooperstock, in Proceedings of the 14th International Conference 
on General Relativity and Gravitation, Florence, Italy, Aug 1995. 

G. P. Perry and F.I. Cooperstock, in Proceedings of the 6th Canadian Conference on 
General Relativity and Relativistic Astrophysics, Fredericton, Canada, May 1996. 
V.S. Manko and N.R. Sibgatullin, Phys. Lett. A 168, 343 (1992); J. Math. Phys. 34, 
170 (1993). 

V.S. Manko and N.R. Sibgatullin, Class. Quantum Grav. 10, 1390 (1993). 
W. Kinnersley, J. Math. Phys. 14, 651 (1973). 
W. Simon, J. Math. Phys. 25, 1035 (1984). 

C. Hoenselaers and Z. Perjes, Class. Quantum Grav. 7, 1819 (1990). 

J. L. Synge, Relativity: The General Theory, (North Holland, Amsterdam, 1966), 
pp. 313. 

In cases A-C, Eq. ( |5.2|) has been solved to a precision of \K\ < 10~ 50 using highly 
refined values of z?,- 

It should be noted that having identical roots a± = 0:2 is not sufficient for identifying 
critically charged bodies even in the Weyl-class. The Curzon particle is such an object 
with «i = a?2 but it is not necessarily critically charged. See Refs. |29| and 0]. 

H. E. J. Curzon, Proc. Lond. Math. Soc. 23, 477 (1924). 



18 



Aguirregabiria J M, Chamorro A, Suinaga J and Vishveshwara C V 1995 Class. Q 
turn Grav. 12 699 



19 



FIGURES 



FIG. 1. Schematic of two Reissner-Nordstrom black holes in cylindrical coordinates. The thick 
lines are the Weyl "rods" which show the locations of the event horizon surfaces. 

FIG. 2. Schematic of a Reissner-Nordstrom black hole and a Reissner-Nordstrom superextreme 
body. The dotted line is a "complex Weyl rod". The intersection of the "rod" with the z-axis is 
defined as the coordinate position of body 2. 
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TABLES 

TABLE I. The values of the parameters in the parameterizations of papers I and II is shown 
given the Weyl-Class values. 



Weyl-Class Paper I Paper II 



mi = 8 


mi = 14.17 


mi = 3.58 


9i =8 


9i = 14.17 


9i = 3.58 


m 2 = 3 


rfi2 = —3.17 


m 2 = 7.42 


92 = 3 


92 = -3.17 


92 = 7.42 


zi = -7 


zi = -2.02 


h = -4.7 


z 2 = 7 


z 2 = 2.02 


z 2 = 4.7 
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Body 1 




Body 2 



Figure 1: Schematic of two Reissner-Nordstrom black holes in cylindrical 
coordinates. The thick lines are the Weyl "rods" which show the locations 
of the event horizon surfaces. 
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o- 



z 2 

i_ 



Body 2 ^ Im (a 3 ) = — Im (a 4 ) 



Z 2 = —Re (a 3 ) = —Re (a 4 ) 



Figure 2: Schematic of a Reissner-Nordstrom black hole and a Reissner- 
Nordstrom superextreme body. The dotted line is a "complex Weyl rod". 
The intersection of the "rod" with the z-axis is defined as the coordinate 
position of body 2. 
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